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1 Introduction 



Future high precision measurements of anisotropics of CMB are potentially very pow- 
erful tools for cosmological studies. They may present accurate data on the values of 
the basic cosmological parameters, as well as information about structure formation 
at an early linear stage. There is a very rich literature on this subject; for a review 
and a partial list of references see e.g. 111 addition to angular variations of 

CMB temperature, primordial density fluctuations induce also, as a secondary effect, 
a polarization of the radiation. A good introduction to the subject can be found in the 



pioneering papers 0-||, while a more up-to-date development is given in refs. [10 



18fl . Our approach here is slightly different from that in the last papers and, though 



the results mostly coincide or agree, there are some essential differences which we will 



specify in Sees. 2 and 3 (see also our earlier paper |T9[| ). 

Measurements of the polarization of CMB will permit to obtain an additional 
information which may help to resolve ambiguities in extracting cosmological param- 
eters from the observational data (for a recent discussion of these ambiguities see 
e.g. ref. p0| ). In particular, polarization is quite sensitive to the presence of tensor 
perturbations (gravitational waves) and, as was shown in refs. [[K|-[|14[] a deviation 
from zero of the so called pseudoscalar or "magnetic" part of polarization would be 
unambiguous indication of the presence of gravitational waves. In what follows we 



will further elaborate this issue and generalize some of the results of ref. [16|. 

Geometrical properties of polarization field were actively studied in the recent 
years |fr3|-||TB[|. Important features, which characterize geometry and topology of 
polarization, are the types of singularities of the vector flux lines tangent to the 
direction of maximum polarization |19| (for an earlier paper see ref. [^T|). This vector 
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is parallel to one of the eigenvectors of the Stokes matrix and have the magnitude equal 
to the magnitude of polarization. The points where this vector vanishes, so that no 
direction is determined, are singular points on the polarization map and the character 
of the singularity to a large extend determines the relief of the polarization map, even 
rather away from the points of vanishing polarization. The maps simulated in different 
papers, including ours (see Fig. 4), demonstrate that the behavior (topology) of the 
flux lines of the polarization field at the points where polarization is non- vanishing, is 
determined by the type of the singularity at vanishing polarization. Hence one may 
make a conclusion about the types of the singular points studying polarization maps 
in the regions where polarization is measurable. The properties of these singular 
points are discussed in detail in Sees. 4 and 5. 

Statistical properties of the anisotropy of CMB temperature and CMB polariza- 
tion are of primary importance for an understanding of their origin. At the present day 
only two mechanisms of generation of primordial density perturbations are known: 
inflation and topological defects. Simplest inflationary scenarios predict Gaussian 
perturbations that results in the Gaussianity of the CMB temperature fluctuations 
and polarization at the surface of last scattering. Tests of a Gaussian nature of the 
temperature fluctuations of CMB, together with a study of its polarization, are im- 
portant probes of inflation. A possible method of testing a Gaussianity of CMB is 
a study of statistics of the singular points discussed in Sec. 4.4. The same statisti- 
cal tests may also be useful for discrimination of signal from noise because Gaussian 
distributions of both temperature fluctuation and polarization generated by noise are 
rather improbable. 

Statistics of global geometrical properties of CMB maps, both for temperature 
fluctuations and polarization, can be efficiently studied with Minkowski Functionals 
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(MF) [^2], ^3[. Applications of this approach to the temperature fluctuations and to 
polarization has been discussed in refs. ||24||-||2q| and refs. |21, [19|], respectively. In a 
recent paper |^] genus statistic of simulated polarization maps was analyzed. In Sec. 
6 we apply all three MF's, which are known for two dimensional maps, for the study 
of statistical properties of CMB polarization. 

This paper is organized as follow: The basic concepts and notations are introduced 
in Sec. 2. In Sec. 3 a local description of polarization is compared with a nonlocal 
one. In Sees. 4 and 5 the classification of singular points is discussed and various 
approaches to this problem are compared. The properties of Minkowski Functionals 
are considered in Sec. 6. We finish with Sec. 7 where a short discussion of results 
can be found. 



2 Generalities 

Here we will describe some general properties of CMB polarization and present the 
necessary formalism. This section has an introductory character and to some extend 
is described in the literature. We make a simplifying assumption that the relevant 
angular scales are sufficiently small, so that the corresponding part on the sky is 
almost flat. In this approximation the polarization field on the sky can be considered 
as two dimensional field on flat (x, |/)-plane. The photon polarization is described 
by the second rank tensor ay in the plane perpendicular to the photon propagation. 
By definition this tensor is traceless, because the trace part, proportional to unit 
matrix, corresponds to zero polarization and can be absorbed in the total intensity 
of radiation. It is convenient to expand this tensor in terms of the Pauli matrices a a , 
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a = 1,2,3, which form a complete system in 2 x 2 traceless matrix space: 

a = £ a a a (1) 

The parameter £ 2 is equal to the amplitude of circular polarization, which is not 
generated by Thomson scattering (due to parity conservation), so that it is assumed 
usually that ^2 = 0. In this case the matrix a is symmetric and is determined by two 
functions: 

/ 



a = 



(2) 



Q U 
\U -Q 

The functions Q and U depend upon the coordinate frame; they are components of 
the tensor and obey the corresponding tensor transformation law: 

o! l3 = T?Tja kl (3) 

where the coordinate transformation is given by x\ = T^x^. In particular under 
rotation of the coordinate system with 



c s 
-s c 



(4) 



where c = cos0, s = sin0, and is the rotation angle, the parameters Q and U are 
transformed as: 

Q' = Q cos 2(f) + U sin 20 
U' = -Q sin 20 + [/cos 20 (5) 

It is more convenient in many cases to work with invariant quantities or at least with 
vector ones, whose direction on the polarization map is easy to visualize. There are 



the following invariants (or what is the same, scalars) that may be constructed from 
the second rank tensor. First, of course, it is the trace, Tra = an. In the considered 
case it is zero. The second invariant is the determinant of the matrix a, 

det a = Q 2 + U 2 (6) 

The maximum magnitude of polarization is given by VQ 2 + U' 2 . The direction of 
maximum polarization is determined by one of the eigenvectors of the matrix Ojj (see 
e.g. ref. fill). 



These are the well known algebraic invariants which exist in any space dimensions. 
One may construct two more invariants, using vector operator of differentiation. They 
can be chosen as: 

S = didjdij 

P = ekjdkdidij (7) 

where j = 1, 2 and dj = djdxK In terms of Q and U these invariants are expressed 
as: 

S = (d 2 - d 2 2 )Q + 2d 1 2 U 

P= (dl- Q (8) 

The first scalar invariant exists in any space dimension, while the second pseudo-scalar 
one exists only in two dimensional space, because of the presence of antisymmetric 
pseudo-tensor e^j (analogous antisymmetric tensor in higher dimensions D has D 
indices). These quantities S and P coincide, up to a scalar factor, with the introduced 



in refs. |T0], [12| B and E fields. To our opinion it is more natural to denote them as S 



and P to stress their scalar and pseudo-scalar nature and not as electric and magnetic 
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parts of polarization because these quantities have nothing to do with vectors. In this 
sense we agree with the terminology of ref. |28j (see also [13|, [T| 



An important feature of the pseudo-scalar P is that it vanishes if only scalar 
perturbations induces polarization in CMB. In this case the Stokes matrix can be 
written in terms of derivatives of one scalar function: 



a 



ij 



{2didj - 6 tJ d k d k ) * (9) 



It is straightforward to check that indeed P = 0. We do not share the opinion 
and/or terminology of refs. [HI 117 where it is stated that the corresponding field 



does not possess a curl. As has been shown in ref. |T9[ this is not true and generically 



the eigenvectors of the Stokes matrix are not curless. The validity of this general 
statement can be verified on simple examples. It means in particular that the flux lines 
of the direction of maximum polarization may have a nonzero vorticity in contrast to 
the statement of refs. |14|, [17| . 



If tensor perturbations are non-vanishing, the polarization matrix has the general 
form determined by two independent functions. As is well known, an arbitrary three 
dimensional vector can be expanded in terms of scalar and vector potentials as 

V = grad $ + curl A (10) 

In two dimensions an arbitrary vector can be expressed as derivatives of a scalar and 
a pseudoscalar: 

Vj = + e jk d k $ 2 (11) 

In direct analogy to that, an arbitrary traceless symmetric 2 x 2- matrix can be pre- 
sented in terms of scalar and pseudo-scalar potentials as: 

ay = (2didj - <%<9 2 ) + {e lk 8 k d j + e jk d k di) $ (12) 
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Of course now the pseudo-scalar P defined in eq. (0) does not vanish and this property 
permits to observe possible tensor perturbations by measurement of CMB polariza- 
tion pi|- ||16|| . If = then the scalar S vanishes. Unfortunately it would not mean 
that tensor perturbations dominate because they contribute both into $ and ^. 



3 Local and nonlocal description of polarization 

It is an interesting observational problem which quantity is easier to measure in a 
noisy background, a differential local or an integrated global one. As was stated 
in ref. [|15j a measurement of an integrated quantity would be much more robust. 



And correspondingly the field variables S and P (or in notations of papers [1^1 -[12 



E and B) were expressed as integrals over all or a part of the sky. We think that 
the answer to the question on the best observational strategy very much depends on 
the properties of the noise. For example if the noise in polarization field of CMB 
is created by point-like sources, chaotically distributed on the sky with the mean 
separation larger than the resolution of the antenna, then the measurement of local 
differential quantities, as e.g. direct measurements of S and P given by eq. (^) seems 
easier. However there may be sources of the noise that would be easier to suppress if 
one measures a quantity which is averaged over whole (which is not possible) or a part 
of the sky. To this end we re-derive the expressions for integrated S and P (or E and 
B) presented in ref. [[U|. The derivation is presented in great detail because of some 



disagreement with ref. ||16|| . The results are very close but we show that the window 
function may have a more general form even for the same choice of normalization 
function N(l 2 ), defined below. 
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Let us first define the Fourier transformed fields: 

Q(f)= Jd 2 ye- ir *Q(y) (13) 

and the similar one for U. The Fourier transformed scalar and pseudo-scalar fields 
can be written as 

S~ N (l) = N(l 2 ) J d 2 ye-^\Q(y) cos 20; + U(y) sin 20,] 

P N (f) = N(l 2 ) J d 2 ye-^[U(y) cos 20, - Q(y) sin 20,] (14) 

where 0, is the polar angle in the plane of Fourier coordinates I. The scalar function 
N(l 2 ) is arbitrary. It preserves scalar or pseudo-scalar property of S and P. For 
the definition (|7|) one has to choose N(l 2 ) = I 2 . The definition used in ref. |16| is 
N(l 2 ) = 1. This means that a non-locality is introduced in the coordinate space by 
the inverse Laplace operator, 1/d 2 , that is by the Green's function of the Laplacian. 

Now we can make the inverse Fourier transform to obtain the functions Sn and 
Pjy in coordinate space: 

Sn(x) = I J^N(l 2 ) J dV r(aH0 IQ(V) cos 20, + U{y) sin 20,] 
Pn(x) = J j^y 2 N(l 2 ) J d 2 ye lf ^ \-Q(y) sin20, + U(y) cos20,] (15) 

where 0, is the angle between the vector I and some fixed direction; it is convenient 
to choose the latter as the direction of vector x, so that 0, = 0^,. 

Integration over directions of vector I can be done explicitly To simplify the 
notations let us introduce 

p = x-y (16) 
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and three angles (pi p , <p px , and <p x i between the directions of the indicated vectors. 
Evidently 



i P + 4>px + 4>xi 



o 



(17) 



The angular integral is reduced to 



277 



d(j)i p e 



Up cos <f>i p 



(A cos 20/p + B sin 20, 



i pi 



The second term 



(19) 



where the coefficient functions A and B do not depend on cj) 
vanishes, while the first one gives 

d(l) lp e Upcos ^ cos20 ip = -2tt J 2 (Ip) 

Jo 

where J 2 (z) is the Bessel function (see e.g. p9[| ). 

The integration over magnitude of / depends upon the form of the function N(l 2 
and the result is a function of the magnitude of the vector p: 



dllN{l 2 )J 2 {lp) = F N {p) 



(20) 



For the particular case of N(l 2 ) = 1 chosen in ref. [|T(j] the integral can be taken as 
follows. It is formally divergent so some regularization procedure should be applied. 
This can be achieved by introducing a small imaginary part to / to ensure convergence 
(in other words, we have to shift the contour of integration to the upper /-half-plane). 
Using the relation |29| : 



zJ 2 (z) = Ji(z) - zJ[(z) 



(21) 



and integrating by parts, we obtain: 



F 1 (p) = - / dzzJ 2 {z) = — 2 / dzMz) 
p Jo p z L Jo 



- zJAz) 



(22) 
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Now taking all the contributions together we obtain: 

^ poo p2n 

Sn(x) — — / dppF N (p) / d(f) [Q (x — p) cos 2(J) + U (x — p) sin 20] 
Z7T Jo jo 

I /"GO /■27T 

P JV (f ) = — / dppF N (p) / d0 [-Q (x - p) sin 20 + U (x - p) cos 20] (23) 

27T JO JO 

For the particular case of F N (p) = Fi (p) = 2/p 2 considered in ref. [T(J we obtain 



almost the same result as the quoted paper with the only difference that we do not 
see any reason to assume that the window function F\(p) = 2/p 2 should be taken zero 
at p = 0. Anyhow, this difference has zero measure and does not have any impact 
on the value of the integrals (f23"|). Hence it may be disregarded. What, as we think, 



is more essential is the statement of ref. [161 that in order to avoid difficult for even 



impossible) integration of the data over the whole sky one may use a modified window 
function: 

2 

? 

with the function g(p) subject to the condition: 



F sz (p) = -g{p) + ^ f dp'p'g(p') (24) 
p z JO 



dppg(p) = (25) 

where the last integral is taken over all the sky. 

We believe that any window function can be used and no additional conditions 
are necessary. To show that we calculate the functions S^(x) and Pn{x) for the 
particular case of scalar perturbations when the Stokes matrix is given by expres- 
sion ([5J). Calculation of derivatives in polar coordinates is straightforward and after 
some algebra we obtain: 



I poo pZn I 

Sn(x) = 7^J q dppW(p) J o d(f) Nf p>p (x - p) 
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where sub-p or sub-0 means differentiation with respect to the corresponding variable 
and W(p) is an arbitrary window function. 

One can see from the second of these expressions that indeed P vanishes for any 
window function. Thus to prove the absence of tensor perturbations, one should 
either observe vanishing of the local quantity P(x) given by eq. ([8p or of the nonlocal 
one given by eq. ( [26]) with an arbitrary convenient window function W(p). Which 
method would be more efficient depends upon the properties of the noise. 



4 Singular points in polarization maps 
4.1 Introductory remarks 

Polarization state of photons of CMB can be described by the direction of maximum 
polarization and its magnitude; the former is parallel to the eigen- vector of the Stokes 



matrix and the latter is equal to \/Q 2 + U 2 . Polarization maps simulated in different 
papers present the corresponding vector field on two-dimensional plane. For the 
analysis of the flux lines on this map it is very important to know the properties of 
the singular points of this vector field. It has been done in ref . |19| (see also |2lJ ) . 



Another approach was taken in refs |15|, [16[] where the properties of the flux lines 



were analyzed in terms of basis functions of tensor spherical harmonics (see figs. 1 
in refs. |15|, |16[). However the flux lines of these basic functions are quite different 
from the behavior of the flux lines of the polarization vector. Of course the analysis 
in terms of tensor harmonics and the behavior of Stokes parameters Q and U can be 
used for description of polarization maps but the analysis in terms of eigenvectors of 
the Stokes matrix permits to make a more direct description of the properties of the 
polarization field. Possible types of different singular points as well as their statistical 
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distributions may bring a new piece of information about properties of CMB. Of 
course a measurement of CMB polarization near the point where it vanishes is a very 
difficult observational problem. However it is not necessary to go exactly to the point 
where Q 2 + U 2 = 0. The type of the singularity can be determined by the pattern 
created by the flux lines in the region where polarization is non-vanishing (see an 
example of simulated polarization map below in Fig. 4). 

The analysis of singular points of the polarization vector field was performed in 
ref. [19|| and f2l| . It was found in [19 that their types do not fit the well known 



classification of singular points of vector fields in the standard theory of dynamical 
systems. Due to non-analytic behavior of the eigen-vectors near the the zero points, 
Q 2 + U 2 = 0, the separatrices end at the singularity, while in the usual case they 
smoothly continue through these points. These unusual behavior, found in our pa- 
per [T9|], is well observed in the polarization maps simulated in refs. pi and in the 



maps of our paper below(Figs. 4,5). In this and the next section we present a further 



development of the analysis of our previous shorter paper [19 



4.2 Basic equations 

The eigenvectors of polarization matrix (|2|) are: 

n + ~ {U,X-Q} 
n~ ~ {-U,\ + Q} (27) 

where A = y/Q 2 + U 2 is the magnitude of the eigenvalue and the vectors n ± corre- 
spond respectively to the positive and negative eigenvalues, ±A. The vector n + is 
parallel to the direction of the maximum polarization, while n ~ goes along the direc- 
tion of the minimal polarization. This is evident in the basis of eigenvectors where 
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the polarization matrix is diagonal, a = diag{A, —A}. The total intensity of light 
polarized along n ^ is given by I± = Iq ± A. Thus the intensity along n + is bigger. 

For defmiteness we will consider the field of directions of the vector n + and the 
singular points in this field. The problem of singular points of a vector field V is 
investigated for the case when the direction of this two-dimensional vector field with 
the components [x(t),y(t)} is governed by the equation 

dy = F 1 (x,y) 
dx F 2 (x,y) 

Singularities may appear if simultaneously both functions F 12 (a;,?/) vanish. In this 
case the conditions of uniqueness of the solution of the differential equation is not 
fulfilled and more than one integral curve may pass through the same point. The 
standard theory is developed for the case when the functions F 12 are analytic near 
these zeroes, and their first order Taylor expansion has the form: 

Fj = aj(x - x ) + bj(y - y ) (29) 

The following three singular points are possible in this case: knots, saddles, and foci 



(see e.g. POfl). The separatrices of the solutions are two intersecting lines, which are 
simply straight lines in the linear approximation. 

However in the case of polarization vector field the basic equation has the form: 



dy_ _n^__ X-Q 
dx nt U 



(30) 



The singular points may appear, as above, if both numerator and denominator vanish. 
It is equivalent to the condition Q = U = 0. An essential difference to the standard 
case is that now the numerator is not analytic near zero. This fact results in a quite 
different behavior of the integral curves near these points. The standard theory is not 
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applicable to this case and below we will investigate the structure of solutions in the 
vicinity of these points directly. We assume that the functions Q and U are analytic 
near the points where Q = U = 0, so that they can be expanded as: 

Q q±x + q 2 y 

U ~ u\x + u 2 y (31) 
For the sake of brevity we assume that Q and U vanish at x = y = 0. 

4.3 Types of singular points 

It is convenient to introduce the new coordinates: 

f = qix + q 2 y 

1] = mx + u 2 y (32) 

Since this coordinate transformation corresponds to a rotation and rescaling of the 
coordinates, the forms of singular points would remain the same. Now we introduce 
polar coordinates on the plane (£, rf): 

£ = r cos <f), rj = r sin <ft (33) 

In these coordinates the equation (|30|) is rewritten as 

dlnr = N^ q 2 t 3 + fa - 2u 2 )t 2 - (q 2 + 2 Ul )t - q x 
d(j) D~ u 2 t 3 + (m + 2q 2 )t 2 + (2 qi - u 2 )t - u x [ ' 

where t = tan(0/2). 

In the general case the denominator D has three roots tj, j = 1, 2, 3. Without loss 
of generality we may assume that u 2 = 1. Then these roots satisfy the conditions: 
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ht 2 + t 2 t 3 + t 3 tx = 2gi - 1, 
h+t 2 + t 3 = -(u 1 + 2q 2 ) (35) 

The integration of equation (B3j) becomes straightforward if we expand the r.h.s. in 
elementary fractions: 

dcj) j f ~ t i 

where, as one can easily see, Bj = N(tj)/ (tj — t k ){tj — t/), none of j, k, I are equal to 
any of the others. It is straightforward to verify that 

B _ (l+t 2 t 3 )(l + tf) 2 
1 2(h - t 2 )(h - t 3 ) 1 ' 

Remaining parameters B 2 and B 3 are obtained by cyclic permutations. 

Since d In r/d<f) = (dlnr /dt)(l + t 2 )/2 the equation can be finally rewritten as 

and the integration becomes straightforward. The corresponding solution is: 

3 



r = tq 



(l+^nc*-*^ (39) 



where tq is an arbitrary constant and the powers Uj are 



"> = r+5 (40) 

with the constants Bj given by eq. fl3~7|) . It can be checked that Vj satisfy the following 
conditions: 

S>i = -i> ( 41 ) 

j 
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3 3 



E^ = 4(E^ + IK] =91, (42) 



2 , 

3 \ 3 3 



" 8(tx - t 2 f{t 2 - t 3 ) 2 (t 3 - Uf (1 + tlt2)(1 + hh){1 + hh) (43) 



The last three factors in equation ([43) are proportional to the determinant d = 
q x u 2 - q 2 ui. 

(1 + t x t 2 ) (1 + t 2 t 3 ) (1 + * 3 ti) = 2(giu 2 - Mig 2 )/^ = 1d/u 2 2 (44) 

If all the roots tj are real, then the sign of the product YG Uj is the same as the sign of 
the determinant d. If however one of the roots, e.g. t\, is real and the other two are 
complex conjugate, the sign of the determinant and of the product fl43|) are opposite. 

Now we can make the classification of the singular points. Let us first consider 
the case when all the roots tj are real. The behavior of the solution is determined 
by the signs of the powers Uj. Due to equation ( fill) at least one of the powers Uj 
must be negative. To see what other signs are possible let us assume (without loss of 
generality) that 

h > h > t 3 (45) 

In this case the following sign relations are valid: 

sign[^] = sign [-(1 + t 2 t 3 )) , 
sign[z/ 2 ] = sign [(1 + ttt 3 )] , 

sign[z/ 3 ] = sign [-(1+ Ma)] (46) 

If e.g. t 3 > 0, the following signs of Vj are realized (— , +, — ). If t 3 < but t 2 > 0, 
then v 3 < and one or both v\ and v 2 are negative. They cannot both be positive 
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because if (1 + tits) > 0, then (1 +M3) > too and v\ < 0. Analogously, in the case 
t\ > and t2 < the set of signs (— , +, +) for any sequence of Uj is impossible. In 
the case when all tj are negative, the sign pattern is (— , +, — ). Thus only two sign 
combinations for Uj are possible: (—,—,—) and (—,—,+). The first one is realized 
when d < in accordance with expressions ( f43|) and ([0]). If the determinant is 
positive, then the signs of Uj are (—,—,+). 

In the case when d < the solution does not pass through zero in the vicinity of 
the singular point. Its behavior is similar to the usual saddle with the only difference 
that there are three and not four, as in the usual case, linear asymptotes/separatrices 
(see Fig. la). We will also call it a "saddle". The fact that in our case separatrices 
are not continued through the singular point, in contrast to the usual singularities 
is related to the non-analytic behavior of the equation fl5D| ) due to the square root 
singularity. 

If d > 0, then the sign pattern is (— , — , +) and the solution vanishes along one 
of the directions and tends to infinity along the other two. The form of the solution 
is quite different from the standard ones. The field line cannot be continued along 
<p = <pi into <p = <pi + it as can be done in the usual case. We will call this type of 
singularity a "beak" (see Fig. lb). 

If only one of the roots tj is real and the other two are complex conjugate, the 
solution has the form: 

- = {f + 1) I t - t 2 | 4Re " 2 exp (4/3 Im v 2 ) (t - h) 2ui (47) 
^0 

where (3 = tan _1 [Imt 2 /(t — Re^)]- The real root v\ is negative, as is seen from eq. 
([43|) and thus r does not vanish in vicinity of such singular point. The flux lines of 
the polarization field for this case are presented in Fig. lc. This type of singularity 
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a b c 

Figure 1: Flux lines for three different types of singular points: (a) saddle, (b) beak, 
and (c) comet. Dashed lines show peculiar solutions (separatrices) . 

can be called a "comet". This case is realized when determinant d is positive. 

4.4 Probabilities of various types of singularities 

The relative weights of different singular points was calculated in the following way. 
It is evident that the probability of "saddles" is 50% because saddles appears if and 
only if d < 0. The probability of comets and beaks was found numerically from the 
conditions that d > and there is a single real root of the equation D = (for comets) 
or there are three real roots (for beaks), where D is the denominator in the expression 
(P^). The probability of appearance of saddles, beaks and comets for random choice 
of qi, 92, Ui, and «2, is correspondingly W s = 0.500, Wt, ~ 0.116, W c ~ 0.384. 

One can also estimate the number density of the singular points in the following 
way (see e.g. || |21|). All singular points correspond to the case when both Q = 
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and U = 0. The number density of these points is proportional to 

dQdU = \d\dxdy (48) 

and thus the density is given by the average value of the determinate, d = qiU2 — q2Ui- 
It can be shown that saddles make 50% of all singular points (n s ) = 0.5(n), where 
n is the number density of all singular points. Calculations of the number density of 
beaks and comets are more complicated and should be done numerically. According 
to our estimates the surface densities for beaks and comets are correspondingly (n&) ~ 
0.052(n) and (n c ) ~ 0.448(n). Deviations from these and above found numbers for 
W s ,b, c m ay signal deviations from Gaussian nature of perturbations. 



5 Comparison of two classification methods 

As we have mentioned above the first classification of singular points in CMB polar- 
ization map was performed in ref. plf , where the equation 

^ = 9. (49) 
dx U 1 } 

was used to describe the behavior of the flux lines of the "vector" V = [U, Q]. However 

under coordinate transformations U and Q are not transformed as components of a 

vector but as components of a second rank tensor in accordance with eq. (|5|). Because 

of that the maps of the flux lines of the " vector" V would not be invariant with respect 

to rotation of the coordinate system. In a fixed reference frame there could be three 



possible types of singular points in accordance with the standard classification [pO 
knots, foci, and saddles. This is so because both functions Q and U generically should 
be analytic near the point where both of them vanish. The types of singularities in 
these classification scheme depend in particular upon the sign of the determinant 
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d = q\U2 — q2U\ introduced above. If d > 0, then the singular points are saddles, 
while for d < there may be both foci and knots. 

It can be shown that d is invariant under coordinate rotation. Thus saddles retain 
their identity in different coordinate frame. On the other hand, foci and knots may 
transform into each other under rotation. So the topology of the map of flux lines 
of the "vector" V would look differently in different coordinate systems, though the 
positions of the singular points evidently remains the same. Moreover the flux lines 
of another possible "vector", W = [Q,U] are in a sense complementary to those of 
V. The relevant determinant changes sign and thus saddles of V correspond to foci 
and knots of W and vice versa. 

In contrast to classification based on V or W, the description of the polarization 
map in terms of eigenvectors of the Stokes matrix, considered in the previous sections 
and in ref . [ I9| , is invariant with respect to coordinate transformations, so that types 
of the singular points do not change under rotation. As we have argued in the 
previous section, positive d gives rise to either to beaks or comets, while negative 
d could produce only saddles. Thus a saddle in the flux lines of corresponds 
to a knot or a focus on \^-map and to a saddle on H^-map. A comet or a beak on 
^(+)-map both correspond to a saddle on V-map or to either a knot or a focus on 
W-m&p. The latter relations may be different in different coordinate frames because, 
as we mentioned above, knots and foci transforms into each other under rotation, 
while beaks and comets remain the same (this is also true for saddles on n^ + ^-map). 
These statements are illustrated on Figs. 1-3 where examples of the maps for n^ + ' 
and for V are presented. 

The simulated 500*500 pixels, 5° x 5° degree map of the CMB polarization field for 
the Standard CDM model with the HWFM resolution 0.3° is presented in Fig. 4. Each 
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— * 

Figure 2: Transformation of the flux lines of "vector" V due to rotation of coordinates 
near the saddle type singularity of n^ + \ plotted in Fig. 1(a). The set of maps from left 
to right and from the top to the bottom corresponds to rotations with respect to the 
first one by the angles: = 5°, 10°, 15°, 20°, 22.5°, 25.6°, 30°, 34°, 49°, 57°, and 70°. 
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Figure 3: Transformation of the flux lines of the "vector" V near the beak (upper 
panel) and comet (lower panel) type singularities of f^ + \ plotted in Fig. 1(b) and 1(c) 
respectively. The maps from left to right correspond to rotation by = 45° and 90° 
with respect to the first one. 
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Figure 4: Simulated map of CMB polarization vector field n^ + ' (see Sec. 5 for details). 
Solid lines show the flux line behavior near singular points where polarization vanishes. 

vector represents orientation of the linear polarization with the angle ^tan _1 (f//Q) 
counterclockwise to the positive direction of the X axes. Length of each vector is 
proportional to \fQ 2 + U 2 . For visual clarity we use only 50*50 pixels. Solid lines 
represent the behavior of the field in the vicinity of non-polarized points. In this 
picture we have 4 saddles, 3 comets and 1 beak. For comparison of two approaches 
described in this section the same random Gaussian realization as for Fig. 4 but for 
the pseudo vector V is presented in Fig. 5. As we discussed above the types of the 
singularities in this case are different. There are 4 saddles which correspond to comets 
or beaks in the polarization map of Fig. 4 as well as 2 knots and 2 foci. The latter 
can only be the saddles in the polarization pattern of Fig. 4. 
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Figure 5: Simulated map of the flux lines of "vector" V (see sec. 5 for details). Solid 
lines show the flux line behaviour near singular points where polarization vanishes. 
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6 Minkowski Functionals. 



In this section we will consider an application of Minkowski Functionals (MF) for 
description of statistical properties of CMB polarization. For two dimensional maps 
the following three MF's are of interest: area (A), length (L), and genus (G). As usu- 
ally, we will consider the magnitude of polarization, ^/Q' 1 + U' 2 as a two-dimensional 
surface in a three-dimensional space. If we cut this surface at different levels P t , then 
the area of the map will be divided into two parts: the area where polarization is 
above the threshold P t and the area where P < P t . Three functions characterize this 
division: 

1. "Area", A, is the fraction of the area of the map where P > P t \ 

2. "Length", L, is the surface density of length of the boundary between the 
fractions with P > P t and P < P t ; 

3. "Genus", G, is the number density of isolated highly polarized regions minus 
the number density of isolated weakly polarized regions (equivalent to the Euler 
characteristic). 

All three Minkowski Functionals, A, L, and G, are evidently functions of Pt- 

Minkowski Functionals (MF) have some mathematical properties that make them 
especially convenient for our task. They are translationally and rotationally invari- 
ant, additive, and have simple and intuitive geometrical meaning. Moreover, it 



was shown |2^, |23| that morphological properties of any pattern in D-dimensional 
space, which possess the above mentioned characteristics, can be fully described by 
D + 1 Minkowski Functionals. The CMB polarization, that we consider, is a two- 
dimensional random field, so the introduced above three MF's would fully charac- 

26 



terize these properties. MF's are additive with respect to isolated parts of the sky 
that makes them potentially very powerful for patchy coverage. These properties 
make MF's a convenient tool for checking the polarization patterns for a presence of 
a non-Gaussian component. 

For statistical analysis of CMB polarization field it is convenient to express the 
Stokes parameters Q and U in terms of scalar and pseudo-scalar potentials \I/ and $ 
as given by eq. (|T^): 

Q = $ u - $ 22 - 2*12, U = 2$i2 + *22 - *22 (50) 

where §>ij, ^ are the second derivatives of $ and * in a Cartesian coordinate system 
in the small (flat) part of the sky, as is described above. The potentials $ and * can 
be considered as two independent two-dimensional Gaussian fields. For description 
of their statistical properties it is convenient to use their Fourier representation: 



2tt J 



= _L f C^{k)e~ l ^dk (51) 
2tt J 

Statistical properties of the fields Q, U and their derivatives are characterized by 
dispersions: 



n* 1 
n* 1 



rZn roo pZtt roo 

/ d6 cos 2 (26) / k 5 \C$(k)\ 2 dk + / d6sm 2 (29) / k 5 \Cl(k)\ 2 dk 
Jo Jo Jo Jo 

dO cos 2 (2#) / k b \C^{k)\ 2 dk + / d6sm 2 {26) / k 5 \C$(k)\ 2 dk 
Jo Jo Jo 



Oq = Ou = ^K* + ^o,*], ( 52 ) 

rOO roo 

<j*= k 5 \Ct,(k)\ 2 dk, a 2 ^= k 5 \C*(k)\ 2 dk, 
Jo ' Jo 



a 



oo roo 
2 _ / lJ\n n.\\2jU J* — I 1.7 in fL\\2 
1,* 



k 7 \C^(k)\ 2 dk, a 2 = k 7 \C^(k)\ 2 dk, 
ii ' Jo 

OO _ _ /"OO 

~dk. 



roo roo 

a 2 2q> = k 9 \C^(k)\ 2 dk, a% = k 9 \Cy(k)\ 
Jo ' Jo 



27 



Expressions for the correlators of first and second derivatives, namely, Qi, Qij, Ui, Uij, 
through a\ $ , o\^, of $ and o\^, are presented below. 

Two important characteristics of the polarization are the autocorrelation radius, 
r c , and the coefficient of cross correlation between functions Q and U and their second 
derivatives, 7: 

rl = §44^ (53) 



7 = + a ^ (54) 



It is convenient to use instead of Q and U and their derivatives dimensionless 
variables, normalized to their dispersions: 

Q U 

Q = r. u = 



(4* + 4*) /2 J (4* + 4*) /2 



^* / rr rr \ 

Ui = (55) 



(*?,* + *?,.) /2 ^(4* + 4*)/ 2 

) ; Uij 1 

\J (4* + 4*) /2 y (4* + 4*) /2 

Below we calculate Minkowski Functionals using these dimensionless variables. It is 
noteworthy that while Qi = dQ/dx\ its dimensionless analogue, is not equal to 
dq/dx 1 but differs from it by the constant factor r c . Thus, for example, 

qi = rAq, = r%djq (56) 

6.1 Area 

The simplest Minkowski Functional, that is the fraction of area of the map with 
p > p t , can be calculated using the joint probability distribution function (PDF) 
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B(q,u) for values of q and u only. Since q and u are Gaussian fields, this function 
can be written using correlators for q and u and their dispersions: 



< q >=< u >=< qu >= 0, < q 2 >=< u 2 >= 1 (57) 

and we have 

1 / q 2 u 2 \ 
B(q,u)dqdu = — exp I -— - — j dqdu , (58) 

The PDF of the amplitude p = y/q 2 + u 2 is: 

B(p)dp = pexp (-y J dp (59) 



and the relative area of the map where p > p t is: 



00 / 2 \ 

A( Pt )= jB(p)dp = expl-^\ 



6.2 Length 

The condition 



(60) 



p(x,y)=p t (61) 



describes a curve in two-dimensional plane. The mean density of the length of this 
curve can be found from the expressions 

dl 



dp(x, y) = (Vp dl) = ^p 2 x + p 2 y -, (62) 
L(Pt) = ^P 2 x +P 2 y )B(p t ) (63) 

where {^Jp 2 + ) must be found for a given p = p t . 
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To do this we should use the PDF B(p t ,p x , p y ). For Gaussian q and u this function 
would also have a Gaussian form and to find it we need correlators of q and u and 
their first derivatives: 



< qqi >=< uu,i >=< qui >=< uq t >= 0, 

< $ >=< m >=< q l u j >= 0, (64) 

I, 

< QiQj >=< UiUj >= -dij 
where is the Kronecker symbol. Using this correlators we obtain 

B(p,p x ,p y )dpdp x dpy = -pexp (-p 2 /2 - p 2 x - p 2 ) dpdp x dp y , (65) 



and, finally, 



(Vi* + i*> = T- < 66 ) 



6.3 Genus 

As is well known, Genus, G{p t ) is expressed through the conditional mean value of 
determinant det(pij) = P11P22 — v\i under conditions p > p t , Pi = 0: 

= ^J dp J U d Pij det (Pij) B (P^Pi = °>P«) = —I^ J d P B (P)( A (p) > (68) 

T ' c Pt -oo ^ T ' c pt 

where B(p,pi = 0,Pij) is the conditional PDF of its arguments and (A(p) ) is the 
conditional mean of det(pij). To find (A(p) ) we should know the conditional mean 
value and dispersions of p^ which can be expressed through the conditionless mean 
value and non- vanishing correlators of qij and My, namely: 

7 

< qqu >=< uua >= —^Sij, 
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2 1 7 + 5p 

< Qu > 



< lit, > 



16 1 + p 2 ' 
1 5 + 7p 2 
16 1 + p 2 ' 



2 1 1 + 3// 
< gii?22 >=< <?12 >= 77 



16 1 + p 2 ' 



2 1 3 + p 2 

< Mll^22 > = < U 12 >= — 



16 1 + p 2 ' 



< <7n«i2 >=< u u q 12 > 
< qnun >=< u 22 q 12 >= 



1 1 



16 1 + p 2 ' 
1 -1+p 2 



16 1 + p 2 ' 

where p 2 = cr^^/o- 2 ^. These correlators are linked to correlators pij through 
expressions 



Q 



u 



1 



Pij = -qij + -uij + Wij, Wij = -{<i,<ij + 
p p p 

The required determinant (A(p) ) can be written as 



U,U 



l LA J 



PiPj) 



2 2 
(A(p) ) = — (qnq 2 2 ~ q\l ) + ^-(Wll«22 ~ «12 ) + -r(?ll«22 + Wll<?22 ) 

f p *p *p 

+ -(Wlig 2 2 + 911^22 } + -(U>11«22 + U U W 22 ) 
p p 

and for conditional mean values we have 



(911922 - <?i 2 2 > = ^(? 2 - 1) 

7 2 

M2-"? 2 ) = -j(g 2 - 1) 

7 2 

(<?ii«22 + Uuq 22 ) = 2— qu 

i 2 q 



(w u q 22 + 9iiW 22 ) 
(wuu 22 + M11W22) 



4 p 

7 2 -u 
4 p 
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Finally 



(A(p)) = ^(p 2 -3), 
G(Pt) = ^rM ~ ^ ( 73 ) 



More cumbersome method for the calculation of the Genus was used in ref. |21 | 



where the calculations have been done under incorrect assumptions about the corre- 
lations of Pij. In spite of that, the final result obtained there remains true. 

All three MF's have very simple analytical forms and can be used, together with 
other methods, to test the Gaussianity of CMB. 



As was noted in ref. []3T| and later discussed in refs. and |19J the condition 
G(p pr ) = defines the critical amplitude p pr for which the regions of high polarization 
percolate. According to eq. (|73|) this critical amplitude of polarization is 



V = Ppr 



that was tested in simulations made in ref. ]2T |. 

As has been mentioned above, G(p t ) is equal to the number of isolated regions 
above p t minus the number of isolated regions below p t . For p t C 1 the percolation 
takes place, i ~ 1 and almost all map is occupied by the one region with p > p t . 
This means that only small isolated spots around the points p = can be separated 
and the mean number density of such points, N Q , is 

N = 1 - G(0) = (74) 



As we can see from eqs. (|53|) , (54) and (56) MF's have an universal form for any 
relation between the potentials $ and Therefore, Minkowski Functionals for the 
CMB polarization field do not depend on the type of perturbations and are determined 
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by the statistics of the signal only. To discriminate the scalar and pseudo-scalar 
perturbations and to estimate the parameter p given by eq. (69), a more detailed 
investigation of statistics of Qij,Uij and pij needs to be performed. This implies, 
for example, direct estimates of some of the correlators (69) which would be a very 
difficult problem because of a finite resolution of the future polarization maps. 

7 Summary and discussion 

The main goal of this paper is to present a critical discussion of the known methods of 
analysis of statistical and geometrical properties of CMB polarization and to propose 
some new additional methods. 

We have compared local and nonlocal descriptions of the polarization with dif- 
ferent expressions for the field functionals S and P (see Sec. 3) in the case of noisy 
data. We conclude that the optimal observational strategy strongly depends upon 
the properties of the noise. For a sufficiently rare sources of noise on the sky, a mea- 
surement of local quantities S and P (0) which are obtained by differentiation of the 
Stokes parameters seems to be more favorable. However an isotropic noise with zero 
average may be easier to suppress by taking an average value over a part of the sky. 

We established a classification of the singular points of the flux lines of the eigen- 
vector r?( + ) of the polarization (Stokes) matrix. We have found that there are three 
possible types of singularities at the points where polarization vanishes: saddles, 
comets, and beaks. They are different from the singularities known in the standard 
theory of two dimensional vector fields. The former, as is well known, are saddles, 
knots, and foci (even though the name "saddle" is the same in both cases the behavior 
of the flux lines are quite different). The reason for this difference is a square root 
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singularity of the vector at the point where polarization vanishes. We have 
compared our method which was earlier presented in ref. |19| and described in detail 



here, with other methods [21, 15, 16 



It is shown that in the case of Gaussian primordial fluctuations (as predicted by 
inflationary cosmology) the saddles make 50% of all singular points, beaks make 5.2%, 
and comets make 44.8%. Deviations from these numbers may signal deviations from 
a Gaussian nature of perturbations or may indicate an existence of a non-Gaussian 
noise in observational data. We realize that a measurement of polarization near the 
points where it vanishes is a formidable observational problem. But the patterns on 
polarization maps created by the flux lines of n^ + ' in the regions where polarization 
is non-zero is to a large extent determined by the types of nearest singularities. Thus 
one may hope that measurements of polarization on the patches where it is sufficiently 
large would permit to determine the singularity types. 

Finally we applied Minkowski Functionals to the analysis of the statistic of the 
geometrical quantities on the polarization maps. It is shown that it provides a sensitive 
test of Gaussianity of CMB. 
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